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^ i We have mapped the C I ^Pi -^^Pq line at 492 GHz in three 

Qv^ I molecular clouds immersed in weak ultraviolet radiation fields, TMC-1, 

0> I L134N, and IC 5146. In all three clouds, the C I peak T^ ~ 1 K, 

1-5 ' with very small dispersion. The spatial C I distribution is extended 

and rather smooth. The J = 2— >1 transitions of CO isotopomers were 
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^^ ■ observed at the same angular resolution as C I. The C I peak TJ is 

c/3 . typically a third of the peak T4 of ^^CO J = 2— >1, and the C I emission 

is usually more extended than emission in ^^CO or C^^^O J=2— >1. The 
C I linewidth is close to the ^^CO J = 2^1 linewidth, larger than 
X ' the C^^O J = 2^1 hne width, and smaller than the ^^CO J = 2^1 

C^ ' linewidth. The shapes of these lines occasionally differ significantly, 

probably because of the combined effects of differing opacities and 
the physical separation of the line forming regions. The uniformity of 
the C I peak T^ is remarkable for a line in the Wien portion of the 
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Planck function and indicates a very uniform excitation temperature. 
This uniformity is best explained if the line if opaque and thermalized. 
If so, the CI line probes kinetic temperature in clouds exposed to 
low ultraviolet fluxes. This conclusion has significant implications for 
the thermal balance in such clouds. At Ay ~2, these clouds have a 
remarkably constant temperature from place to place and from cloud to 
cloud (7.9 lb 0.8 K). Photodissociation region models of clouds immersed 
in the mean interstellar radiation field tend to predict stronger lines 
than we see, but this may be an artifact of assumptions about the 
temperature. 

Subject headings: ISM: abundances — ISM: atoms — ISM: individual 
(TMC-1, L134N, L183, IC 5146) — ISM: molecules 



1. INTRODUCTION 

Neutral atomic carbon (C°) plays an important role in interstellar physics and 
chemistry. C° was recognized long ago as a potentially important coolant for the 
interstellar medium (Werner 1970; Penston 1970). Furthermore, the abundance of 
C° may reflect the physical structure of molecular clouds (e.g. Stutzki et al. 1988). 

Most observations of the submillimeter lines of atomic carbon (C I) have been 
carried out toward Hll regions or reflection nebulae, i.e., toward giant molecular 
clouds with nearby O and B stars. Such regions, where the ultraviolet (UV) 
radiation is strong, are good targets to test PDR ("photodissociation region" or 
"photon-dominated region") models. Recent PDR models can explain the amount 
of C° relative to CO rather successfully (van Dishoeck & Black 1988; Hollenbach, 
Takahashi, & Tielens 1991). PDR models for uniform clouds predict that C° exists 
in the surface layer sandwiched between layers where C+ and CO are the dominant 
carriers of gas-phase carbon (see Keene et al. 1997 for a review). Observationally, 
however, the C I distribution is often similar to the CO and ^^CO distribution. 
Most steady-state models predict a very low abundance for neutral carbon in the 
interior of the cloud (C°/CO<0.04 for Ay >2.5-4, Jansen et al. 1995, Sternberg 
and Dalgarno 1995), but Le Bourlot et al. (1993) have argued that a second 
solution to the chemical equilibrium equations exists at low densities in which the 



neutral carbon abundance is substantially higher (C°/CO=0. 1-0.3). Because of 
strong emission from the PDR, these models cannot be tested in dense, strongly 
irradiated molecular cores. In principle, one might test them in clouds with low UV 
fluxes. If substantial C° exists in the interiors, it may contribute to the formation 
of carbon chain molecules, which are abundant in some dark clouds (Suzuki et al. 
1992). 

Our goals in this work are to determine how C° is distributed in clouds 
exposed to low UV fluxes, and to assess the contribution of such clouds to C I 
emission from galaxies. Molecular clouds which are not illuminated by strong UV 
sources have been observed by several researchers (Keene et al. 1987; Stark & van 
Dishoeck 1994; Ingalls et al. 1994, 1997; Schilke et al. 1995; Stark et al. 1996). 
In these studies, fewer than ten positions toward each were observed with small 
(10-15") beams. The maps of C I emission almost all cover scales of less than 1' 
toward objects with a size of hundreds of square arc minutes. 

In this paper, we present more extensive C I maps and maps of emission 
from various CO isotopomers toward three clouds immersed in weak UV fields. 
TMCl and L134N contain no internal sources of UV radiation. The third source, a 
molecular cloud associated with IC 5146, contains an HII region, but the portion 
of the cloud we observed was 1°-1°5 (20-30 pc) east of the exciting stars. One can 
parametrize the intensity of UV radiation field incident on the clouds in terms of 
the average field in the solar neighborhood, I^/y (Draine 1978). By comparing the 
100 fim intensity observed toward TMCl and L134N to the intensity predicted 
from models of opaque clouds exposed to a UV field with a strength comparable 
to the field in the solar neighborhood (Boulanger & Perault 1988; Bernard et al. 
1990; Hollenbach et al. 1991), we infer that Ijjv ~ 1 for these two clouds. A similar 
analysis for the region we mapped in IC 5146 leads to the conclusion that If/y ~30. 
By studying these sources with weak incident UV fields, we hope to increase 
our understanding of how UV radiation affects the chemistry and energetics of 
molecular clouds with low column densities or of the outer portions of more opaque 
molecular structures. 



2. OBSERVATIONS 

Observations were carried out with the 10.4 m telescope of the Caltech 
Submilhmeter Observatory (CSO)^ between 1994 January and 1994 April, and in 
1994 June, 1995 June, and 1998 March. 

^ The CSO is operated by the Cahfornia Institute of Technology under funding 
from the National Science Foundation, contract AST 96-15025. 

We observed the ground-state fine structure line of neutral carbon ("^Pi ^^Pq, v 
= 492.1607 GHz, Frerking et al. 1989) as well as the J = 2^1 transitions of ^^CO, 
i^CO, C^^O, and C^^O at 219-230 GHz, and the J = 3^2 transitions of ^^CO and 
C^^O at 329 — 346 GHz. For simplicity, we call these two latter frequency bands 
230 and 345 GHz. To make C I and CO observations of large areas on the sky 
in a limited time, we employed a re-imaging device for the Caltech Submillimeter 
Observatory (CSO), with which the 492 GHz receiver illuminates only one panel of 
the main dish (Plume & Jaffe 1995; see also Plume et al. 1994, 1999 for studies of 
molecular clouds with local UV sources with this instrument). With the re- imaging 
device, the measured telescope beam size (6'f,) and Moon efficiency {r\M) are the 
same at both 492 and 230 GHz: Qy, = 2^5 and r]u = 0.8. At 345 GHz, these 
parameters have not been measured but are assumed to be the same as those at 
492 and 230 GHz. Since the emission is very widespread and smooth, riM is the 
appropriate efficiency. Some positions toward L134N were also observed in the 
on-the-fly (OTF) mapping mode (without the re-imaging device), and then the 
data were convolved to match the resolution of the re-imaging device. Because the 
results of the OTF observations are consistent with those of the observations with 
the re-imaging device, the obtained spectra are co-added to improve the S/N ratio. 
For selected positions toward TMC-1, observations without the re-imaging device 
were also carried out on 1994 March 3 and 1994 March 7. For this work, the beam 
size was 15" and the main-beam efficiency was 0.53 at 492 GHz, and the beam size 
was 32" and the main-beam efficiency was 0.76 at 230 GHz. 

A technical description of the 230 GHz, 345 GHz, and 492 GHz front-end 
receivers is given by Kooi et al. (1992), Kooi et al. (1994) and Walker et al. (1992), 
respectively. The SSB system temperature was 700—5000 K, 700—1100 K, and 
180-500 K (SSB) at 492, 345, and 230 GHz, respectively. The receiver backend 
was a 1024 channel acousto-optic spectrometer. The total bandwidth, channel 
width, and spectral resolution were 49.7 MHz, 48.5 kHz, and 140 kHz, respectively. 



C I spectra were binned to 195 kHz (0.12 km s^^) or 390 kHz (0.24 km s^^) 
channels depending on the S/N ratio. Data were obtained by position switching 
against a position shown to have no emission. The spacing interval employed for 
mapping was 3'. The line temperature scale was calibrated by using the standard 
chopper-wheel method. Throughout this paper, we give intensities on the corrected 
antenna temperature (T^) scale (see Kutner & Ulich 1981 for definition). The 
calibration assumes that the signal and image sidebands have equal gain. 

We observed the classical "dark" clouds, TMC-1 and L134N (a.k.a L183) and 
a portion of the cloud associated with IC 5146. In addition, we use data toward 
regions with local UV sources, Orion A (this work), Cep A, NGC 2024, S140, and 
W3 (Plume et al. 1999), for comparison. The reference center, off position, and 
observed lines for each object are listed in Table 1. 



3. RESULTS 

For all the observations, we removed a first order baseline, determined an 
integrated intensity from the area under the line, and fit a Gaussian to determine 
peak T4. While some lines were non-Gaussian, the fits always provided a good 
estimate of the peak T4. 



3.1. TMC-1 

The distance to TMC-1 is 140 pc (Ehas 1978; Cernicharo et al. 1985). At this 
distance, the beamsize of the re-imager (2f5) corresponds to 0.10 pc. Figures 1(a) 
and (b) show the distribution of the C I and ^^CO J = 2— >1 spectra on the sky. 
The map center (Table 1) is the cyanopolyyne peak (Churchwell, Winnewisser, & 
Walmsley 1978), and the map size is 18' (EW) x 33' (NS) (or 0.7 pc x 1.3 pc). 
The peak line temperature decreases slowly away from the map center, except to 
the south. The mean JT^(C I)dV is 1.7ib0.3 K km s^^. (Throughout the paper, 
the uncertainty refiects the standard deviation of the distribution, not the error 
in the mean.) The corresponding values for T^ are 1.1 it 0.3 K (Fig. 12). The 
"TMC-1 ridge", which is a ~ 2'-wide NW-SE ridge passing through the map center 
in high-density tracers such as NH3, CCS and HC3N (Little et al. 1979; Hirahara 



et al. 1992), is not prominent in C I or ^^CO. Figure 2 shows the distribution of the 
jTXdV along (a) Aa = 0" and (b) A6 = 0" . The distribution of /r|(C I)dV, like 
that of ^^CO J = 2^1, is fairly flat across the source, while the C^^O J = 2^1 
and C^^O J = 1^0 (Langer et al. 1995) peak more sharply at the map center. 

The C I and ^^CO J T^dV correlate very poorly in this source. From all the 
points in the re-imager map of TMC-1, we obtained a least-squares fit to a relation 
between the / TXdV of: 

jTXiC I)dV= (1.1±0.4) + (0.19 ± 0.11) / Tl{^^CO) dV 

where the correlation coefficient was 0.32. 

This poor correlation contrasts with the strong correlations (correlation 
coefficient 0.85) found in large-scale observations of UV-illuminated molecular 
clouds (Plume et al. 1999) and the even stronger correlations in the cloud cores 
themselves (Keene et al. 1997). In our large-scale TMCl map, the C I profile does 
not always resemble the ^^CO profile. To see the difference in the line shape more 
clearly, we averaged the spectra of three positions: (Aa, A6) = (0", —1260"), (0", 
-1080"), and (0", -900"). Figure 3 shows the average spectra of CO J = 2^1, 
i3co J = 2^1, C^^O J = 2^1, and C 1. The C 1 profile is flat-topped with almost 
constant temperature from 5.3 to 7 km s~^. The ^^CO J=2^1 profile has a sharp 
peak at 5.5 km s~^, a broad plateau, and then a second peak at ~6.8 km s^^. 

If some of the C I emission arises in high column density molecular cores and 
if the C° abundance depends on the local chemistry, changes in that chemistry 
might affect the strength of the C I emission. To investigate this possibility, we 
measured the ^ T^^dV of C I and ^^CO at the cyanopolyyne peak (0", 0") and at 
the ammonia peak (Act, A5) = (-290", 380") (e.g., Hirahara et al. 1992; Pratap 
et al. 1997). Because the ammonia peak position does not match our observation 
grid, we average the spectra at (-360", 360") and at (-180", 360"). S T\{C I)dV 
is 2.0 and 1.8 K km s~^ at the cyanopolyyne and ammonia peaks, respectively, 
while the ^^CO ^ T\dV is 4.6 and 4.9 K km s~^. There is no significant peak of 
the C I line at the cyanopolyyne peak; if deep C° contributes to long chain carbon 
molecules, it is not visible in this line. 

To see the C I line temperature and spectral variation on smaller scales, we 
observed the detailed line emission distribution across the TMC-1 ridge with the 
fully illuminated CSO (beamsize 15" for C I, 30" for C^^O 2^1). Figure 4 shows 



C I and C^^O J = 2^1 spectra. The grid spacing was \/2 x 20" from NE to SW. 
The C I and C^^O jT^dV distribution is rather featureless on this (1') scale. The 
mean value of fT^{C T)dV is 2.18 K km s~^ with a standard deviation of ±0.36 
K km s-i. For C^^O J=2^1, we find jTXdV = 1.31 ± 0.19 K km s^^ The C I 
profile is broader than that of C^^O (AV(C I)/AV(C^^O))=1.9±0.5. Our results 
are consistent with those of Schilke et al. (1995), who observed the region with the 
same telescope but with a somewhat different reference center. 



3.2. L134N 

The distance to L134N is llOiblO pc (Franco 1989). The cloud is isolated, and 
it shows no star forming activity in far-infrared observations (Sargent et al. 1983; 
IRAS Point Source Catalog 1985). L134N has been extensively observed in various 
molecular lines (Swade 1989). Phillips & Huggins (1981) previously detected C I at 
one position and Keene (1995) has made a 20' long strip-map of C I in this source. 
Figure 5 shows spectra obtained in an area of 12' (EW) x 30' (NS) (0.4 pc x 1.0 
pc). The beamsize using the re-imager corresponds to 0.08 pc. For C I, the jT^dV 
distribution is remarkably uniform; the mean and standard deviation are 1.5 it 0.5 
K km s-^ The mean T^ is 1.1 ± 0.2 K (Fig. 12). 

Figures 6(a) and (b) compare the C I and CO isotopomer spectra at the 
map center. The C I profile peaks at 6.8 km s~^ and has a prominent blue wing 
extending to Ylsr= km s-\ The ^^CO, C^^O, and C^^O J = 2^1 and ^^CO 
and C^^O J=3^2 lines differ significantly in their widths, shapes, and central 
velocities. These differences imply that line opacity and/or isotopic fractionation 
effects influence the shapes of most of these lines. 

Figure 7 shows the jT^dV distribution along Aa = 180". The C I distribution 
is very flat, varying by only 20% along a strip more than half a degree (1 pc) long. 
As in TMCl, the correlation between the C I and ^^CO jT^dV, based on the 
entire map, is poor: 

JTX{C I)dV = (0.94±0.29) + (0.14 ± 0.07) / TH^^CO) dV 

where the correlation coefficient was 0.41. 

Stark et al. (1996) observed six positions along a lO'-long East- West strip in 



L134N with the James Clerk Maxwell Telescope. They found that the distribution 
of JTX{C I)dV is similar to those of ^^CO J = 3^2 and J = 2^1 but not to 
C^^O J = 2^1. Our result shows that the distribution of ^ T\{C I)dV is flatter 
than those of ^^CO and C^^O J = 2^1. Ratios of ^^CO/C I J TXdV vary from 1.8 
to 3.0 along this strip. The distribution of j T^{C I)dV is also substantially more 
extended than the emission in the many molecular transitions observed by Swade 
(1989), including C^^O J =1^0, CS J = 2^1, NH3 J,K = I, 1, Hi3C0+ J = 
1^0, and C3H2 Jk-k+ = lio ^loi- 



3.3. Cases with Local Ultraviolet Sources 

3.3.1. IC 5146 

IC 5146 (distance = 1.0 kpc, Walker 1959) is an open cluster with an H 11 
region. At 1 kpc distance, the beamsize of the re-imager corresponds to 0.7 pc. A 
neutral cloud associated with the H 11 regions has been observed in molecular lines 
by Dobashi et al. (1992), Lada et al. (1994), and Kramer et al. (1999). Figure 8 
shows spectra observed in an area of 39' (EW) x 21' (NS) (11.8 pc x 6.4 pc). The 
reference center is the position of IC 5146. We have observed a region ~ 60' — 90' 
away from the H 11 region. Observations of this portion of IC 5146 serve as a bridge 
between the study of clouds immersed in the weak ambient UV field in the outer 
Galaxy and regions illuminated by prominent local UV sources. JT^{C I)dV is 
2.0ib0.9 K km s~^ over the mapped region. The peak T^ is typically about 1 K 
(0.9 lb 0.3, Fig. 12). Despite our estimates of somewhat higher Ifjv in this region, 
the lines are not significantly stronger. 

Figure 9 shows the J TXdV distribution along A6 = 1380". The C 1 and ^^CO 
intensities are poorly correlated. From all the IC 5146 data, we obtain 

JTX{C I)dV= (1.18±0.53) + (0.22 ±0.12) / T^i^CO) dV 

with a correlation coefficient of 0.44. 



3.3.2. Other Sources 

In order to understand the differences between C I emission from sources with 
different luv, we make use of large-scale C I and CO isotopomer maps of high-/;7V' 
sources. The data include points from the maps of Cep A, NGC 2024, S140, and 
W3 made with the re-imager by Plume et al. (1999). In addition, as part of the 
present work, we have mapped a 15' x 25' region in Orion A with the re-imaging 
device. Figure 10 shows the J T^dV maps of Orion A in the C I ^Pi -^^Pq and 
^^CO J = 2^1 emission. luv for regions 2-15 arcmin from 9^ C Orion ranges 
from 10 to 10^ (Stacey et al. 1993). The T| is much higher and more variable than 
that in the clouds with low UV fluxes: mean TJ is 3.8 it 1.1. White & Sandell 
(1995) observed a smaller part of the Orion A cloud in C I with the James Clerk 
Maxwell Telescope. With its smaller beam, they found that the distribution of C I 
is often different from that of ^^CO. Nonetheless, Keene et al. (1997) show that 
C I emission correlates well with ^^CO emission in Orion at moderate intensity 
levels. In our data, the C I emission on the map is doubly peaked around (Aa, A6) 
= (0", —300") and (0", 0"), while the ^^CO emission does not show double peaks. 
Our map is consistent with the map of White & Sandell (1995), if the difference in 
beamsize is taken into account. 

Figure 11 plots the C I jT^dV versus the ^^CO jT^dV for the three dark 
clouds (low-Iuv) as well as for the five high-Ijjv molecular clouds. A least-squares 
fit to the data from all of the high-Ijjv sources gives: 

JTX{C I)dV= (3.8±0.3) + (0.25 ± 0.02) / TX{^^CO) dV 

with a correlation coefficient of 0.75, clearly a stronger correlation than found for 
the low-/(7y clouds. The dark cloud points cluster among the points from the 
high-J^y clouds with the lowest C I and ^^CO J = 2^1 ^ Tj^dV . The overlapping 
points from the high-Jf/y clouds mostly represent regions near the edges of the 
clouds and often those parts of the clouds with the smallest amounts of incident 
UV radiation (Plume et al. 1999). Even in aggregate, the dark clouds show little 
or no correlation of C I and ^^CO j T^^dV . 
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3.4. Intensity and Linewidth 

Since the mean T^ is similar in all the low-uv clouds, including IC5146, we 
combined them to find the distribution in Fig. 12; the mean r| is 1.0 it 0.3. This is 
a far lower value and a tighter distribution than is seen in clouds exposed to higher 
UV fluxes. 

In contrast, the linewidths, relative to CO isotopes, are similar to those in 
high-Iuv clouds. Plume et al. (1999) obtained (AV(C l)/AY{'^^CO)) = 1.08±0.3 
and (AV(C-'^^0)/AV(^^CO)) = 0.78ib0.2 for massive star forming regions with 
local UV. For TMCl, L134N, and IC 5146, we measured the linewidth ratio in the 
same way to be (AV(C I)/AV(i3C0)) = 0.97±0.21 and (AV(Ci80)/AV(i3CO)) = 
0.58ib0.16, which are similar to the results for sources with stronger local UV. Here, 
we used spectra toward (Act, A6) = (0", -720")- (0", 180") in TMC-1, (0", 0") 
and (180", -540")-(180", 360") in L134N, and (-5580", 1380") - (-5220", 1380") 
and (-4320", 1380") - (-3960", 1380") in IC 5146 where both lines were observed 
and the line shape is not so much different from a Gaussian profile. The absolute 
values of the linewidth are AV(C I) = 1.3±0.4 km s-^ AV(CiSO) = 0.7±0.2 km 
s~\ and AV(i3cO) = 1.3±0.2 km s"! for TMC-1 and L134N. Those for IC 5146 
are AV(C I) = 1.5±0.4 km s-\ AV(CiSO) = 1.3±0.5 km s-^, and AV(i3co) = 
1.8ib0.4 km s^^. When we compare the C I linewidth with the ^^CO linewidth, we 
obtain (AV(C I)/AV(12qo 2^1)) = 0.62±0.05 for (Aa, A6) = (0", -720")-(0", 
180") in TMC-1. 



4. DISCUSSION 

We have detected widespread, remarkably uniform C I emission from TMCl 
and L134N, and from a portion of the IC 5146 cloud far from the OB star cluster. 
The C I emission is more extended than the emission in the J = 2^1 transitions of 
^^CO and C^^O and the strength of the C I line correlates poorly with the jT^dV 
for these two CO isotopomers. In this section, we examine the implications of these 
observations of clouds with low-Ijjv for the origin of the C I emission, for cloud 
chemistry, and for the thermal structure of dark clouds. 
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4.1. Photodissociation Regions Explain the Observed C I Emission 

In our previous paper (Plume et al. 1999), we studied the C 1 '^Pi ^^Pq 
emission from giant molecular clouds bathed in radiation fields 10-1000 times 
stronger than the mean interstellar field in the solar neighborhood. We showed 
that this emission arises from a neutral atomic carbon layer in a photodissociation 
region at the cloud surface. The measured column densities and intensities agree 
well with the values predicted by the theoretical models. Both the independence of 
C I jT^dV from the strength of the UV field and the higher C I/^^CO ratios at 
the cloud edges also support the conclusion that the C I emission comes from the 
photodissociated material at the atomic/molecular interface. We would now like 
to use the C 1 and CO isotopomer maps of TMCl and L134N, clouds with luv ~ 
1, to examine the case for PDR's as the source for the C I emission from clouds 
immersed in lower UV fields than those incident on the CMC's. 

Figure 11 plots the C I J T^dV versus the strength of the ^^CO J=2^1 fine 
for the positions we have observed in TMCl and L134N and, for comparison, the 
positions in IC 5146 and other clouds illuminated by stronger local UV fields. The 
figure shows that the C I jT^dV in TMCl and L134N are comparable to those 
observed in IC 5146, and comparable to those at the weakest points observed in 
the extended molecular clouds near OB stars observed by Plume et al. (1999). In 
this figure, we also present theoretical C I and ^^CO J T^dV calculated on the 
basis of column densities predicted by models of translucent clouds illuminated by 
UV radiation (van Dishoeck and Black 1988). The models include the effects of 
isotope-selective photodissociation and isotopic fractionation and assume a carbon 
depletion factor, 6c of 0.4. The open squares represent models T3 to T6, clouds 
illuminated by luv = 1 and total cloud extinction. Ay*, from 2.0 to 5.1. The open 
circles represent models 14 to 17, clouds illuminated by I[7y=10 and total cloud 
extinction. Ay*, from 2.7 to 9.0. We calculated the J T^dV from the model C° and 
^^CO column densities using a Large Velocity Gradient (LVC) (Scoville & Solomon 
1974; Goldreich & Kwan 1974) non-LTE excitation and radiative transfer code 
to account, in a crude way, for excitation and opacity effects. We used the rate 
coefficients for collisional excitation of C I from Schroder et al. (1991). The C° 
and ^^CO column densities in the models do not depend strongly on the kinetic 
temperature assumed by van Dishoeck and Black (see their models T6A and 
T6B). Uncertainties in predicted intensities from our LVG calculations therefore 
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depend directly on the assumed kinetic temperature rather than indirectly through 
linewidth-dependent photochemical effects. In this rough comparison, we have used 
Tx=10 K to be consistent with measures of core temperatures in dark clouds (e.g. 
Benson & Myers 1980). We set the column density per unit velocity in the radiative 
transfer model equal to the C° or ^^CO column density from the theoretical model 
divided by the typical linewidth in the clouds. The results in Figure 11, where 
we have multiplied the model jT^dV by the beam efficiency to compare to the 
observations, show that the van Dishoeck and Black models for I[/y=l correctly 
predict the C I J T^dV we observe in the maps of TMCl and L134N. The observed 

tot 

V • 



^^CO jT^dV are comparable to the values for the models with the largest A*"* 



The morphology of the C I and molecular emission from the clouds with I(/y ~ 
1 also supports the idea that the C I line arises in material at the cloud surface. In 
both L134N and TMCl, the C I ! T\dV distribution is broader and flatter than 
that of ^^CO and significantly broader than the C^^O J=2^1 distribution. 



4.2. Why Don't C I and Molecular Line Intensities Correlate in 

'Lq-w-Ijjv Clouds? 

If the C I emission arises in a thin layer near the cloud surface and emission 
from molecular lines less opaque than the lowest few transitions of ^^CO arise 
throughout the bulk of the cloud, there is no a priori reason to expect the C I 
and ^^CO or C^^O J T^dV to correlate. Nevertheless, the C I and CO isotopomer 
/ T^dV correlate astonishingly well in high column density cloud cores illuminated 
by powerful UV fields like M17 SW (Keene et al. 1997). There is a clear difference 
in the jT^dV correlation between high-Ijjy clouds (correlation coefficient = 0.85 
for C I vs. i^co jTXdV Plume et al. 1999) and the low-Iuv dark clouds studied 
here (correlation coefficient = 0.32 and 0.41 for C I vs. ^^CO in TMCl and L134N, 
respectively). 

For cloud cores. Plume et al. (1999) point out the similarity of the 
clump /inter clump structure needed to produce the observed extent of the C II 158 
/im line emission to the kind of structure that would produce the high degree of 
spatial correlation between C I and ^^CO seen by Keene et al. (1997); in both cases 
the / T^dV refiect the number of clumps in the beam. Plume et al. argue that in 
extended CMC material, unlike the high column density cores, the C I lines arise 
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predominantly in a global layer over the entire face of the cloud. They explain the 
somewhat weaker C I-^^CO jT^dV correlation as a conspiracy of edge effects and 
the correlation of temperature and density variations with total column density. 
The lack of intensity correlations in TMCl and L134N, then, results from an 
absence of temperature variations as well as from higher opacities in C I ^Pi ^^Pq 
and the lowest transitions of -"^'^CO and C-'^^O at low temperature (see below). As 
is so dramatically illustrated in the small-scale C I-molecular line comparisons 
toward the TMCl ridge (Figure 4 and Schilke et al.), the C I observations toward 
clouds illuminated by Ijjv ~1 external fields do not contain any information about 
high column density molecular regions. 



4.3. Thermal Balance: The Meaning of Constant C I Line Brightness 

The most striking feature of the data are the very narrow distributions of T^ 
seen in Fig. 12. Such tight distributions are not seen in any molecular line besides 
CO, which is opaque and thermalized. We suggest that the CI lines are also opaque 
and thermalized; thus the values of T^ provide a probe of T^ in the PDR layer 
with T ~ 1. 

There are both observational and theoretical arguments for a high opacity in 
the C 1 ^Pi ^^Pq line. The large linewidth and, in TMCl, flat-topped shape of 
the C I lines are both indications that the opacity is high. Furthermore, the peak 
intensity does not vary strongly despite significant variations in line width (Figures 
la, 5, and 8). If the lines were not thermalized, T^ would be very sensitive to 
density fluctuations and the tight distribution seen in Fig. 12 would be extremely 
unlikely. To investigate the C 1 opacity more thoroughly, we also computed 
spherical Monte Carlo models (Choi et al. 1995) for photodissociated clouds. The 
temperature and C I abundance versus distance output from a plane parallel PDR 
model with luv=^ and nH+2H2=^0'^ cm~^ was mapped onto a spherical cloud 
(see Pak et al. 1998, van Dishoeck & Black 1987) with a microturbulent velocity 
dispersion of 0.65 km s^^. The model, which has a line center opacity ~2, correctly 
reproduces the linewidths seen in the ^Pi ^^Pq spectra. Assumptions about the 
cloud core temperature (see below) may mean that the C I opacity is even higher. 
Because the critical density for the ^Pi ^^Pq transition is only 800 cm~^ and the 
photon escape probabilities are substantially lower than unity, it is very likely that 
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the transition is thernialized (Tex = Tk)- 

If the C I line is opaque and thermahzed, it provides a probe of Tk deep in 
the PDR layer, where r ~ 1. The distribution of Tk is shown in Fig. 13 for all the 
dark clouds and for Orion A. The distribution for the dark clouds is remarkably 
tight: mean Tk is 7.9 zt 0.8 K, a 10% variation. The high frequency of the C I 
^Pi ^^Pq transition means that it is in the Wien limit for the temperatures we 
derive; thus a tight distribution in T^ (30% variation) becomes even tighter in Tk- 
Small variations in temperature would lead to an enormous fractional change in 
observed T^ (Ti^=8ib3 K implies, for thermahzed gas, a T^ range from 0.17 to 
2.5 K, a factor of 15). In comparison, the range of Tk inferred for Orion is much 
broader, consistent with a lower and more variable opacity; thus the values of Tk 
in Orion are lower limits. 

Do these values of Tk agree with predictions of PDR models for low-J^y? The 
PDR models discussed above produce Tx ~11.5 K in the line forming region. This 
higher Tk results in a ~3 K C I line. The high predicted temperature may be, in 
part, an artifact of the assumption that the asymptotic cloud core temperature was 
10 K or of slight differences between the actual and theoretical thermal balance. 
Recent models by Kaufman et al. (1999) find that Tk drops below 10 K for very 
small ratios of luy to density. Comparison to their predictions of intensity finds 
agreement with our observations only at low luy and quite high density {n > 10^ 
cm~^). Is there other evidence for such low Tk in low-Ifjy clouds? Clemens, Yun, 
&: Heyer (1991) found that 74% of a sample of such clouds have Tk = 8.5 K. 

Why do the C I layers have such a constant temperature? Given the 
inevitability of density fluctuations, what causes the thermal balance to choose 
this temperature? One strong possibility is that neutral carbon itself provides 
the thermostat. The fact that the C I line is in the Wien portion of the Planck 
function at these low Tk means that the cooling power in the 492 GHz line rises 
very rapidly with temperature in the 5-10 K range. In the specific layer where C I 
is present, the higher opacity of the ^^CO line and the low abundance of ^^CO 
may keep the rotational lines of CO isotopomers from contributing strongly to the 
cooling. Indeed models by Hollenbach, Takahashi, & Tielens (1991) indicate that 
C I cooling does dominate in the layer where neutral atomic carbon is abundant. 

If our interpretation is correct, the higher transition of C I, ^P2 —^^Pi at 809 
GHz, will be very weak toward these clouds. An upper limit can be obtained by 
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assuming that the higher transition is also opaque and thermahzed. In this case, 
the radiation temperature (or T^/tjm) is 0.3 K. However, the higher hne is very 
unhkely to be opaque. Even if the levels are thermalized, the optical depth in the 
higher line will be only 0.16 that of the lower line. More realistic estimates can be 
obtained from models of clouds in weak UV fields, combined with LVG calculations. 
Model T6 (n = 10^ cm^^ A'^^ = 5.1 mag) of van Dishoeck & Black (1988) and 
LVG calculation with Tk = 8 K predict the ^Pi ^^Pq and •^P2 ^'^i^i radiation 
temperatures to be 0.9 K and 0.02 K, respectively. Model T6C with higher density 
(n = 10^ cm~^) and LVG calculation with Tk = 8 K predict the ^Pi -^^Pq and 
^i-*2 -^^Pi radiation temperatures to be 1.1 and 0.04 K, respectively. The expected 
ratio of ^P2 ^^Pi to ^Pi ^^Pq is thus less than 0.05. To the authors' knowledge, 
no observations of C I ^P2 -^^Pi toward molecular clouds in weak UV fields have 
been reported in the literature. 



4.4. Overall Cooling by C I 

The C I ^Pi ^^Pq line was expected to play an important role in interstellar 
gas cooling (Penston 1970). The relevant quantity here is -F = J J ludi^d^; because 
F oc h'^jT^dV, C I may be more important than one might guess from the ratios 
of jTXdV. We derive F (C l)/F (^^CO J = 2^1) = 1.9±0.3 for TMC-1 and 
1.0±0.5 for IC 5146. The ratio F (C l)/F (^^CO J = 3^2) in L134N is 0.8±0.2. 
In these dark clouds, the cooling power of the C I ^Pi ^^Pq line is as important as 
low- J ^^CO lines. The flux ratios in the average spectrum of the Galaxy (Wright et 
al. 1991) are F (C l)/F (^^CO J = 2^1) = 2.3 and F (C l)/F (^^CO J = 3^2) 
= 1.6. 



5. Summary 

We have mapped the C I ^Pi ^'^Pq transition in three clouds with low-I[/y 
values at their surfaces. The C I emission is more extended than emission in 
low-J CO isotopomer lines and the intensity is very uniform. C I linewidths are 
significantly larger than C^^O linewidths and comparable to or larger than the 
widths of low-J ^^CO lines. The uniform intensity and large linewidth of the C I 
lines, as well as models of low-I^/y photodissociation regions imply that the C I 
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^Pi ^^Pq transition is thernialized and optically thick. The uniform T^=1.0ib0.3 
K, together with the high frequency of the C I line, means that the kinetic 
temperature in the emitting layer is remarkably constant, T/^(C I layer)=7.9ib0.8 
K. We suggest that the C I emission itself acts as a thermostat in controlling the 
temperature of this layer in low-luv clouds. 

The authors would like to thank Wenbin Li for his help during the observations 
and the CSO staff for their technical assistance. Thanks are also due to Jill 
Knapp, Ken Young, Merce Crosas, and Zeljko Ivezic for letting us use a part 
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AST-9017710 and AST-9530695 to the University of Texas and by the David 
and Lucile Packard Foundation. K. T. was supported by Grants-in-Aid from the 
Ministry of Education, Science, Sports, and Culture of Japan (Nos. 07740180 and 
11440067) and by a grant from the Sumitomo Foundation (No. 950570). 
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Table 1: Reference Center and Off Position 



Source Reference Center 
a(1950) <5(1950) 



Off Position 

Aa A6 Observed Line 



TMC-1 4'^38™38^6 25° 36' 00 



o ocl r\t\ll 



Orion A W"iT''AT -5° 24' 23' 



L134N 15''51™30'' -2° 43' 31" 



// // // 

IC 5146 2l'*51™54'' 47° 01' 12" 



2400" 
1800" 
-660" 



-1800" 
-3000" 



0" C I, 

J = 2^1 of I^CO, I^CO, C^^O 
1800" C I, 

J = 2^1 of I^CO, C^^O 
660" C I, 

J = 2^1 of I^CO, C^^O, C^^O, 
J = 3^2 of C^^O 
0" ^ J = 3^2 of I^CO 
0" C I, 

J = 2^1 of I^CO, i^co, C^^O 



" Off position is in (AAZ, AEL) for this line only 
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Fig. 1.— Line profiles toward TMC-1. The (0,0) position was a(1950)= 4''38™38^6, 
(5(1950)= 25° 36' 00". These data were taken with the re-imager with a 150" beam. 
The intensity scale is T^. (a) C I ^Pi ^'^Pq, rebinned to 0.24 km s^^ channel 
spacing, (b) ^^CO J = 2^1. 

Fig. 2.— Distribution oi J TXdV toward TMC-1, along (a) Aa = 0" (b) A5 = 0". 
Open circles, open boxes, and filled circles represent the jT^dV of ^^CO J = 2^1, 
C^^O J = 2^1, and C I ^Pi -^'^Pq, respectively. The error bars reflect only the rms 
noise level. 

Fig. 3. — Line profiles observed with the re-imaging device toward TMC-1. C I 
spectra were binned to 0.24 km s~^ channels, and some of the CO isotopomer 
spectra were binned to 0.12 km s^^ channels. The spectra are sums of data taken 
at three positions: (Act, A5) = (0", -1260"), (0", -1080"), and (0", -900"). 

Fig. 4.— C I 3Pi ^3Po (left) and C^^O J = 2^1 (right) spectra obtained toward 
TMC-1 with the full CSO telescope (beamsize 15" for C I and 32" for C^^O J=2^1). 

Fig. 5. — C I line profiles toward L134N. These data were taken with the re-imaging 
device with a 150" beam. The intensity scale is T^. The (0,0) position was a(1950)= 
15'*51™30^ (5(1950)= -2° 43' 31" Spectra were binned to 0.24 km s^^. 

Fig. 6. — A comparison of line profiles observed with the re-imaging device at the 
map center for L134N. (a) The J = 2^1 transitions of ""^'^CO, C"*^^, and C^^O, and 
C L C I, C^^O, and C^^O spectra were binned to 0.24 km s^^ channels, (b) The 
J =3^2 transitions of ^^CO and C^^O. The spectra were binned to 0.24 km s^^ 
channels. 

Fig. 7. — Distribution of J T^dV along Aa = 180" in L134N. The symbols are the 
same as those used in Figure 2. 

Fig. 8. — C I line profiles toward IC 5146. The (0,0) position for this map was 
a(1950)= 21^51^54'* (5(1950)= 47° 01' 12". 

Fig. 9.— Distribution of the jTXdV along A6 =1380" in IC 5146. The symbols are 
the same as those used in Figure 2. 
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Fig. 10.— Map of jTXdV of C I ^Pi -^^Pq (solid lines) and I^CO J 2^1 (dotted 
lines) toward Orion A. The contour interval is 3 K km s~^ for C I and 10 K km s~^ 
for I^CO. 



Fig. ll.~ A plot of C I jTXdV against that of ^^CO J=2^1. The filled diamonds, 
triangles, and boxes represent values observed toward TMC-1, L134N, IC 5146, 
respectively. The dots represent points from the CMC's observed for this paper 
and by Plume et al. (1999). The connected open symbols represent C I and ^^CO 
/ T^dV predicted by an LVG model and column densities from the PDR chemistry 
models of van Dishoeck & Black (1988). The open squares represent models T3- 

, tot 



T6, for I[/y=l and total cloud extinction. Ay*, from 2.0 to 5.1. The open circles 



represent models 14-17, for It7V'=10 and total cloud extinction. Ay*, from 2.7 to 9.0. 



Fig. 12. — The lower three panels show the distribution of T| for the three low-Ijjv 
clouds, with the mean and standard deviation indicated. The distribution for all 
three dark clouds is shown in the upper left, while that for Orion (note very different 
scale) is shown in the upper right. 



Fig. 13. — The distribution of Tk, assuming that the lines are opaque and 
thermalized, is shown for all the dark clouds and for Orion A. The assumption 
that the line is opaque is probably not valid in Orion A, so these represent lower 
limits to Tk- 
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